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We investigate a cosmological scenario with three interacting components that includes dark 
matter, dark energy, and radiation in the spatially flat Friedmann-Robertson-Walker universe. We 
introduce a 3-dimensional internal space, the interaction vector Q = (Q x , Qm, Qr) satisfying the 
constraint plane Q x + Q m + Q r = 0, the barotropic index vector 7 = (7:r,7m,7r) and select a 
transversal interaction vector Qt in a sense that Qt ■ 7 = 0. We exactly solve the source equation 
for a linear Qt, that depends on the total energy density and its derivatives up to third order, 
and find all the component energy densities. We obtain a large set of interactions for which the 
source equation admits a power law solution and show its asymptotic stability by constructing the 
Lyapunov function. We apply the x 2 method to the observational Hubble data for constraining the 
cosmic parameters, and analyze the amount of dark energy in the radiation era for the above linear 
Qt. It turns to be that our model fulfills the severe bound of Q x (z ~ 1100) < 0.1 and is consistent 
with the future constraints achievable by Planck and CMBPol experiments. 
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I. INTRODUCTION 



Based on the large and still growing number of astro- 
nomical observations, one can agree that there is an ex- 
otic component called dark energy which represents more 
than the 70% of the total energy of the Universe. Its 
existence has radically changed our standard paradigm 
of cosmology mostly because of its visible effects on the 
current state of the Universe It turns out that dark 
energy is a repulsive fuel characterized by a strong neg- 
ative pressure to overcome the slowing down effect of 
gravity, making the Universe exhibit an accelerated ex- 
pansion state at the present time [l|. This tremendous 
fact has been confirmed by a plethora of observational 
tests such as the high redshift Hubble diagram of type 
la supernovae as standard candles 0] and accurate mea- 
surements of cosmic microwave background anisotropics 
3. According to the current observations, the present- 
day value of the dark energy density is about 120 orders 
of magnitude smaller than the energy scales at the end of 
inflation, so one of the main challenges in the modern cos- 
mology is to understand such deep mismatch. One way 
to alleviate the aforesaid problem is working within the 
context of dynamical dark energy models, leaving aside 
the standard ACDM model. Besides, the necessity of 
a dark matter component comes from astrophysical evi- 
dences of colliding galaxies, gravitational lensing of mass 
distribution or a power spectrum of clustered matter [l[ , 
[H, 0]. The first evidence of dark matter's existence 
stemmed from the studies performed by Zwicky in 1934 
to the Coma cluster of galaxies @ and since its discovery, 
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dark matter has played an essential role for resolving the 
riddle of the missing mass in the Universe. At the present 
moment, the astrophysical observations from the galac- 
tic to the cosmological scales suggest that dark matter 
is a substantial component to the Universe's total mat- 
ter density [B[ and sustain that dark matter represents 
nearly 25% of the total energy matter of the Universe; 
this invisible and nonbaryonic component is the major 
agent responsible for the large-structure formation in the 
Universe [l|, Q|. Motivated to understand more about 
the nature of both dark components, one could consider 
an exchange of energy between themselves, i.e., the dark 
matter not only can feel the presence of the dark en- 
ergy through a gravitational expansion of the Universe 
but also can interact between them 0. A coupling be- 
tween dark energy and dark matter changes the back- 
ground evolution of the dark sector, allowing us to con- 
strain a particular type of interaction and giving rise to 
a richer cosmological dynamics compared with noninter- 
acting models [7j . One way to extend the insight about 
the dark matter-dark energy interacting mechanism is to 
explore a bigger picture in which a third component is 
added, perhaps a weakly interacting radiation term as it 
occurs within a warm inflation paradigm Q. A scenario 
in which dark energy interacts with both dark matter 
and radiation was explored in Q whereas the validity of 
the generalized second law of thermodynamics was stud- 
ied in [l(| without using a particular kind of interaction. 
Other cases correspond to take the third component as 
an unparticle fluid [llj , unparticle fluid in loop quantum 
cosmology |l2|, or a general unspecified fluid As a 
step forward to constraining dark matter and dark energy 
with the physic behind recombination or big-bang nucle- 
osynthesis epochs, a decoupled radiation term was added 
to the interacting dark sector for taking into account the 
stringent bounds related to the behavior of dark energy 



2 



at early times [lij], (l5| . 

Below we develop a model composed of three interact- 
ing fluids and introduce a 3-dimensional internal space 
where the three interaction terms and barotropic indexes 
are viewed as vectors in a vector space. We discuss the 
existence of a transversal interaction and center our find- 
ing in a model with energy exchange proportional to a 
linear combination of the total energy density and its 
derivative up to third order. We also study the stability 
of a power law solution (scaling solution) with the help of 
Lyapunov's theorem. Finally, we perform a cosmological 
constraint using the Hubble data and the severe bounds 
for dark energy at early times. We will use the units 
8nG = 1 and signature (—,+,+,+) for the metric of the 
spacetime. 

II. THE MODEL 

We consider a spatially fiat homogeneous and isotropic 
Universe described by Friedmann-Robertson- Walker 
(FRW) spacetime with a line element given by ds 2 = 
—dt 2 + a 2 (t)(dx 2 + dy 2 + dz 2 ) being a(t) the scale fac- 
tor. The Universe is filled with three interacting fluids, 
namely, dark energy, dark matter, and radiation so that 
the evolution of the FRW Universe is governed by the 
Fricdmann and conservation equations, respectively, 

3H 2 = p = p x + p rn + p r , (1) 

p + 3H(p x +p x + p m + Pm + Pr + Pr) = 0, (2) 

where H = a/a is the Hubble expansion rate. Introduc- 
ing the variable rj = ln(a/ao) 3 , with ao the present value 
of the scale factor and ' = d/drj, Eq. ^ can be recast 

P' = -IxPx - IrnPrn ~ IrPr, (3) 

where 74 = 1 + p^/ pi is the barotropic index of each com- 
ponent with i = {x,m,r} and < j x < j rn < j r . The 
interaction terms ZHQi between the components are in- 
troduced by splitting © into three equations: 

P x +lxPx = Qx- (4) 
P'm + ImPm = Qm, (5) 
Pr + IrPr = Qr, (6) 

where the Qi describe the energy transfer and satisfy the 
condition 

Qx + Qm + Qr=0, (7) 

to recover the whole conservation equation © after hav- 
ing summed Eqs. (d))-®. Also, we assume that all com- 
ponent energy densities are definite positive. 



To investigate the proposed model, we introduce a 3- 
dimensional internal space with an orthonormal vector 
basis {et, e Q , n} defined in the following way: we set the 
coordinate origin at the intersection of the interaction 
plane ([3]) with the vector formed with the barotropic 
indexes 7 = (7^,7 m ,7r), the pair of vectors {e t ,e } 
is contained in the interaction plane, et is orthogonal 
to the vector 7, so 7 • e t = 0, the orthogonal projec- 
tion of the vector 7 on the interaction plane defines 
the direction of the vector e Q and the normal vector 
n = (1,1, l)/\/3 is orthogonal to the interaction plane. 
The interaction terms Qi are also viewed as the compo- 
nents of a vector Q = (Q x ,Qm,Qr) that lives on the 
plane n : Q x + Q m + Q T = 0, meaning that 

Q = qt et + q a e Q , (8) 

where q t and q are the components of the interaction 
vector Q on the plane n and n • Q = 0. 

Taking into account that et is the unique vector of 
the basis with the property of being orthogonal to 7, wc 
adopt this property as a simple criteria for selecting only 
those interactions which are collinear with the aforesaid 
preferred direction in the plane n that we call "transver- 
sal interaction" , so 

Qt = qt e t (9) 

with Qt = qt ensuring that we can always take 7 • Q = 
uniquely. The transversal character of the interaction 
vector © will simplify enough the equations which de- 
termine the component energy densities as we will sec 
below. 




The two basis vectors that span the interaction plane 
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II are given by 



A. Stability analysis 



e t 



(7m - lr, lr ~ lx, lx ~ 7m) 



(10) 



e n = 



(7m + lr - 2"fx, lx +lr - 2j m , "fx + 7m ~ 273) 



V3< 



(11) 



where e 2 = 3 ^ if — 7i] 2 while the vector built with 
the barotropic indexes 7 = 7„n + 7 e lives in the plane 
spanned by the normal n and e Q [see Fig.©]. 

Now, we will construct an interacting three fluid model 
with the transversal interaction ©. After differentiating 
Eq. © and using Eqs. ©-©we have 



P" = llPx + ifnPm + ifpr- 



(12) 



Solving the algebraic system of equations in the 
(Px,Pm,Pr) variables ©, ((3)), and ([T2]). we obtain p x , 
p m , and p r as a function of p, p', and p" only; 

Px = ^ [imlrP + (7m + 7r)p' + p") , (13) 



7x ~ 7r 



[7x7rP+(7, +7rV +/], (14) 



7x - 7n 



[ixlmP + (7 y+p"], (15) 



where A = (i x - 7 m ).(7aj - 7r)-(7m - lr) is the deter- 
minant of that algebraic system of equations. Equations 
(Unj), (HU), and (fT5j) clearly represent the straightforward 
extension of the case studied in Q where an interacting 
two-fluid scenario for the dark sector in the FRW Uni- 
verse was investigated. Following Ref. 0, we replace 
([15]) into ©, (0 into ©, or ([IS]) into © and obtain 
the same third order differential equation, that we call 
"source equation" , for the total energy density; 

p'" + (lx+lm+lr)p" + 
(ixlr + Ixlm + lmlr)p' + (16) 

where its source term Q involves a linear combination of 
the interaction vector components Qi 

Q = IxQmlr + IrQxlrn + ImQrlx- (17) 

By combining the transversal interaction ©, the basis 
vector QJj]) and Eq. {T7]) we find 



S = -A 



(It 



(18) 



Finally, once the transversal interaction Q t is speci- 
fied we obtain the energy density p by solving the source 
equation Eq. (|16p and the component energy densities 
p x , p m , and p r after inserting p into Eqs. (|13p . (1141) . and 



Now, we will investigate the stability of power law 
solutions. To this end, we will assume that the source 
equation ([TB]) admits scaling solutions and subsequently 
we will look for the set of interaction terms that give 
rises to them. The knowledge of power law solutions is 
very useful because it determines the asymptotic behav- 
ior of the effective barotropic index 7 = (ixPx + Impm + 
lrPr)/{Px + Pm + Pr) = —2H/3H 2 , which ranges between 
lx < 7 < lr- These solutions represent a Universe ap- 
proaching to a stationary stage characterized by 7 = 7 S 
and a = £ 2 / 37s ; then the existence of the attractor solu- 
tion 7 S will imply that 7 goes to the asymptotic constant 
value 7 — > is. So on the attractor 



7s 



or 



IxPxs ~t~ IraPras ~t~ IrPrs lx ~t~ ImTmx + 7 r' rx 



Pxs 



J ms ~T~ Prs 



1 



' mx \ 1 rx 



(19) 



{lx - Is) + (lm - ls)r mx + (lr - ls)r r x = 0, (20) 

where we have defined the ratios r mx = p ms /p xs and 
frx = Prs I Pxs- From the positivity of the component en- 
ergy densities, the small value of the ratio r rs = il rs /£l xs , 
the range of the effective barotropic index i x < 7 < j r , 
and tEq. (|20[) . we determine that 7 S ranges between 
lx < Is < lm < lr and the ratios r mx and r rx be- 
come asymptotically constant on the attractor stage, al- 
leviating the cosmic coincidence problem. In the case of 
7 S = 0, we have a final de Sitter regime, H = cte with 

Tmx = -(ix/lm) ~ r rx (lr/lm)- 

To investigate the existence of a power law attractor 
solution a = t 2 / 3js , we use that p' = ~ip, p" = (1 2 — 
l')p, p'" = —(i 3 —3n'+i")p, change the source equation 
([TBI) into 

7" + (lx + Ira + lr - 3jW - V(l) = --, (21) 

P 



V(l) = -(7 - 7s) (7 - 7m) (7 - lr), 



(22) 



and impose both (i) that 7 S be a constant stationary 
solution of Eq. ([2T]l and (ii) the stability condition so 
that 7s is stable. The existence of 7 S implies i' s = 0, 
l'J = and 



Q(ls)=V(is)p, Q(is)<0, 



(23) 



due to Eq. p7]) . 

Taking into account the specific form of the interac- 
tion (|23|) . we will assume separability and the stability 
analysis will be performed for interactions of the form 



Q(p, i, i\ i") = v^JCd, i, i') p 



(24) 



Combining Eqs. (p?T]) and ([2~4"]) . we write the equation 
governing the dynamical evolution of the barotropic in- 
dex in a simpler form 

l" + (lx + lm + lr ~ 37)7' = -V{n)[K - 1], (25) 
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where the function JC fulfills the condition, 

/C( 7 = 7 5 ,7' = 0, 7 " = 0) = l, (26) 

for assuring the existence of the constant stationary so- 
lutions 7 S . Perturbing around the solution 7s by taking 
7 = 7,, + e with |e/7s| <C 1, Eq. (f25j) can be recast as 

e" + (7x+7m+7r-37 s )e' = -P(7 s )[/C 7 e+/Cye'+/C 7 «e"], 

(27) 

where /C 7 , /Cy, and /Cy stand for the partial derivatives 
with respect to 7, 7', and 7", respectively. These deriva- 
tives are evaluated at the point (7 = 7 S , j' s = 0, 7" = 0) 
and we have used that /C( 7 = 7s+e, e') = 1+/C 7 e+/C 7 <e'+ 
/C 7 »e" + C(e 2 , e' 2 , e" 2 ). It turns out that (|27]) can be writ- 
ten as the equation of motion for a dissipative or antidis- 
sipative mechanical system. This resemblance emerges 
from the analogy with the classical potential problem 



d_ 

drj 



-+V(6) 



(28) 



where a = [7 ;c +7 m + 7r-37s + P(7s)£ 7 <]/[l+'P(7s)AC 7 »] 
and the potential 



V(e) - 0- 



(29) 



with = P( 7s )/C 7 /[1 + 7 , (7s)/C 7 »]. In order to assure 
the stability of the scaling solution 7 = 7s , we demand 
that both coefficients a and are positives. So, for any 
transversal interaction Qt leading to a > 0, the poten- 
tial V has a minimum at e = when ft > 0, the func- 
tion inside the square bracket in Eq. (|28j) is a Lyapunov 
function, and the perturbation decreases asymptotically 
reaching e = (attractor) in the limit r\ — > 00; then the 
system is asymptotically stable in the sense of Lyapunov. 
In other words, when the condition (|26[) is fulfilled j s be- 
comes a constant stationary solution of Eq. (|25l) and it 
is stable whenever the stability conditions a > and 
> are satisfied. Note that the interaction (|24p de- 
pends complicatedly on p, p', p" and p'" because the term 
£(7, 7', 7") depends on 7 = -p'/p, 7' = {pi '/ p) 2 - p"/p, 

and 7" = -{?'"/ P) ~ (P'/P) 3 + 3(p'/p)[p'7p - (p'/p) 2 ]; 
then /C contains some nontrivial terms implying that our 
findings are indeed valid even for a broad set of nonlinear 
interactions. 



B. Linear transversal interaction Qt 

Following Q, we assume a transversal interaction Q t , 
with 7 • Qt =0, which is a linear combination of p Xl p c , 
p r , their derivatives up to first order, p. p' . p", and p'", 
so from Eq. (fT5|) we have 

Q = a\p x + a 2 p m + &3Pr + Oiip' 



Using Eqs. (|T3"]). (fT4|) . and (dSJ, we can recast (J3DJ as 

Q = 0lP + 02p'+03p"+04p"'. (31) 

where the new 0j coefficients are written in terms of the 
old cti ones. Eq. (|3"Tj) clearly shows that the most general 
linear interaction only requires a linear combination of p 
and its derivatives up to third order. Combining Eqs. 
and (pTTj) along with p' = —7/9, p" = (7 2 — 7 ')p> an d 
= — (7 3 — 377' + 7")/°, w e obtain that 

ft - 4 ( 7 3 + 7" + 377') + 03 7 2 - 037' - 027 



P 



AC = 



P( 7 ) 



(32) 

Applying the condition (j2"6")l referring to the existence of 
the constant solution 7s to (|32p . we get a constraint 



(33) 



1 -0 47 , 3 +037 2 -027 s =n7 s ), 



for the coefficients 0i, 02, 03, and 04. By solving Eq.([33j) 
for 02 and inserting into Eq. ()31|) . we obtain the final 
form of the effective interaction term 

Q = 0ip + 7 s " 1 [0i-^(7 s )-047 s 3 +037 2 ]p' + 03p" + 04/', 

(34) 

so in the end, the most general linear interaction (|30[) or 
(|34[) only involves four parameters, namely, 0i, 03, 04, 
and 7 S . Replacing (fM]) into Eq. ([2"Tj). we find the poly- 
nomial 7 3 + A7 2 + B7 + C = for constant 7 solutions; 
one of its roots is 7 S , whereas the other two are given by 



A + ~f s 



-1± 4/I 



4C 



7+ 



c 

1-ls 



(35) 



(36) 



where the coefficients A, B, and C turn out to be 

A = J2i H - 03 c = Ixlmjr ~ 01 ^ 



04-1 



04-1 



B = 



Ei^j 7.7, ~ 7s" 1 [0i ~ V{ls) ~ 047 3 + 037s 2 ] 

04-1 



(38) 



The exact solution of the source equation (|T6"|) for the 
general linear transversal interaction (|34l) is given by 



p = b ia - 3 i° + b 2 a- 3 t+ + b z a-^- 



(39) 



whereas the component energy densities are obtained 
from Eqs. Ip jl . ([Tg ]l . and JUJ): 



P, = 7m A 7r [7m7r(&ia' 37s + 6 2 a" 37+ + 6 3 a~ 37 -) 



+ ot^p" + a 7 p"' + «8Pk + a 9 p' m + a w p' r . 



(30) 



-(7m + 7r)(6i7,a" 3Aa + &27+a _37+ + &37-a~ 37 - ) 
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+ (7 S 2 M" 373 + 7 +6 2 a- 37+ + (40) 



Observational constraints on a transversal 
interacting model 



lx - 7r 

A 



[7*7r(M" 373 + b 2 a- 3 ^ + 63a- 37 -) 



-(7* + 7r)(6i7,a" 37a + ^27+a" 37+ + & 3 7-a~ 37 ~) 



+ h 2 s b ia - 3 ^ + 7 +6 2 a- 37 + + 7! 63a- 37 -)], (41) 



Pr = lx A 7m [7m7.(^ia' 373 + fea~ 37+ + ha- 3 ''-) 



'{l x +l m ){b lls a~ 3 ^ +b 2l+ a- 3 ^ +h 1 ^a- 3 ^-) 



+ ( 7 2 M~ 37s + 7+fea- 37 + + 7-&3a~ 37 -)], (42) 

In general, the total energy density in terms of the phys- 
ical quantities such as density parameters is given by 



p = 3H 2 [Ax 3 ^ + Bx 3l+ + Cx 31 - 



(43) 



with x = z + 1 and the z cosmological redshift while the 
integration constants are given by 



A = 



^mo( 7 + - 7m) (7- - 7m) + ^ro(l+ - 7r)(l- ~ 7r) 



(7«-7+)(7»-7-) 



^0(7+ - TxXT- - 7g) 
(7s - 7+) (is - 7-) 



(44) 



B 



^mo(7s - 7m) (7m - 7- ) + ^0r(7s ~ 7r)(7r ~ 7~ ) 

(7s - 7+)(7+ - 7-) 



^xQ^s ~ lx){jx ~ 7-) 

(7s - 7+X7+ - 7-) 



(45) 



Emails - 7m)(7m - 7+) + ^rojls ~ Jr){jr ~ 7+) 

(7s - 7-) (7- - 7+) 



^zo(7s - lx)ilx - 7+) 

(7s - 7-) (7- - 7+) 



(46) 



where fioi = Po«/3i?o are density parameters fulfilling the 
condition Q x q + il r o + O m o = 1 for a spatially flat FRW 
Universe. Additionally, we will choose (7^, 7r „, j r ) — 
(0, 1, 4/3) to recover the three self-conserved cosmic com- 
ponents in the limit of vanishing interaction. 



We will provide a qualitative estimation of the cosmo- 
logical paramaters by constraining them with the Hubble 
data [la]- [l7| and the strict bounds for the behavior of 
dark energy at early times [22| In the former case, 
the statistical analysis is based on the x 2 function of the 
Hubble data which is constructed as (e.g., fl8| |) 



x 2 w = E 



[H(0,z k )-H obs (z k )) 



fe=i 



a(z k y 



(47) 



where stands for cosmological parameters, H b s (z k ) is 
the observational H(z) data at the redshift z k , cr(z k ) is 
the corresponding la uncertainty, and the summation 
is over the 12 observational H(z) data. The Hubble 
function is not integrated over and it is directly related 
with the properties of the dark energy, since its value 
comes from the cosmological observations. Using the ab- 
solute ages of passively evolving galaxies observed at dif- 
ferent redshifts, one obtains the differential ages dz/dt 
and the function H(z) can be measured through the re- 
lation H(z) = -(1 + z)- x dzldt [HI, [Tz|. The data 
H bs(zi) and H b s (z k ) are uncorrelated because they were 
obtained from the observations of galaxies at different 
redshifts. 

From Eq. ([4*3"]) , one finds that the Hubble expansion of 
the model becomes 



H(0\z) = H (Ax 3 ~ ( ° + Bx 4 +Cx : 



(48) 



where A, B, and C are obtained form (jUJ), ([4"5]). and (|4"6"1) . 
respectively. For practical reasons mentioned in the last 
section we have simply selected 7+ = 4/3, 7_ = 1. 

Here, we consider = {Hq, 7s , ft x o, tt m o} plus the con- 
straint on the density parameters to assure the flatness 
condition (fi r o = 1 — f2xo — H m o); then we have four 
independent parameters only. We will take two indepen- 
dent parameters and will fix the other ones along the 
statistic analysis until all parameters have been varied 
and estimated with the x 2 function. Then, for a given 
pair of (01,02), we are going to perform the statistical 
analysis by minimizing the x 2 function to obtain the 
best-fit values of the random variables C = {0\ c ,02c\ 
that correspond to a maximum of Eq. (|47p . More pre- 
cisely, the best-fit parameters C are those values where 
Xmin(@c) leads to the local minimum of the x 2 (0) distri- 
bution. If Xd.o.f. = Xrmn( c)/(N - Tl) < 1 the fit is good 

and the data are consistent with the considered model 
H(z\0). Here, N is the number of data and n is the 
number of parameters (l8j . The variable x 2 is a random 
variable that depends on N and its probability distribu- 
tion is a x 2 distribution for N — n degrees of freedom. 
Besides, 68.3% confidence contours in the 2D plane are 
made of the random data sets that satisfy the inequality 
A X 2 = x 2 {0) - X 2 nm ( e c) < 2.30. The latter equation 
defines a bounded region by a closed area around C in 
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the two-dimensional parameter plane; thus, the la error 
bar can be identified with the distance from the 9 C point 
to the boundary of the two-dimensional parameter plane. 
It can be shown that 95.4% confidence contours with a 
2a error bar in the samples satisfy A^ 2 < 6.17. 




0.15 0.4 0.65 0.4 0.8 1.2 



FIG. 2: Two-dimensional C.L. associated with la,2a for different 6 
planes. 

The two-dimensional C.L. obtained with the standard 
X function for two independent parameters is shown 
in Fig. ©, whereas the estimation of these cosmic 
parameters is briefly summarized in Table (jj). We ob- 
tain that 7 S varies from 10~ 4 to 10~ 3 , so these val- 
ues clearly fulfill the constraint 7 S < 2/3 that assures 
the existence of the accelerated phase of the Universe 
at late times. We find the best fit at (Hq,Q, x o) = 
(74.32km s- x Mpcr x ,0.77) with Xd.o.f = °- 783 b y usin S 
the priors fi m o = 0.2 and j s = 10 -3 . These findings 
show, in broad terms, that the estimated values of H 
and fl X Q are in agreement with the standard ones re- 
ported by the WMAP-7 project The value of n. xa 
is slightly greater than the standard one of 0.7 with a 
discrepancy only of 0.1%. Moreover, we find that us- 
ing the priors (i?o,7e) = (74.20km s" 1 Mpc -1 , 10~ 3 ) the 
best-fit values for the present-day density parameters 
are considerably improved, namely, these in turn give 
(Q x o,Q m0 ) = (0.74,0.23) along with a lower goodness 
condition (xd o f = 0.779). Regarding the estimated val- 
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FIG. 3: Plot of Q x (z) Q m (2), Q r (z), uj e ff(z), and q(z), using the 
best-fit values obtained with the Hubble data for different 8 planes. 
These plots follows the same same order used of in two-dimensional 
C.L. of Fig.l 



2D Confidence level 


Priors 


Best fits 


Xd.o.f 


(O m0 ,7s) = (0.2, 10" 3 ) 


(H ,n x0 ) = (74.32,0.77) 


0.783 


(ii x0 ,is) = (0.70, nn 3 ) 


(tfo.Ono) = (72.28,0.29) 


0.834 


(n x0 ,n m0 ) = (0.7,0.29) 


(ff ,7 s ) = (71.74,10- 4 ) 


0.890 


(#0,7*) = (74.20, 10- 3 ) 


(n x0 ,ft m0 ) = (0.74,0.23) 


0.779 


(H ,n m o) = (74.20,0.24) 


(0^0,7=) = (0-74, 10" 4 ) 


0.855 


(H ,Q x0 ) = (74.20,0.75) 


(7s,ano) = (10- 4 ,0.24) 


0.776 



TABLE I: Observational bounds for the 2D C.L. obtained in Fig. |2j 
by varying two cosmological parameters. 



ues of fi m o, we find that it varies from 0.24 to 0.29, with- 
out showing a significant difference with the standard 
ones [20( | . In performing the statistical analysis, we find 
that Ho G [71.28, 74.32]km s^Mpc -1 so the estimated 
values are met within la C.L. reported by Ricss et al 
[IH, to wit, H = (72.2 ± 3.6)km s" 1 Mpc" 1 . 
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FIG. 4: Plot of logH^O) for z £ [10 2 ,10 5 ] using the best-fit val- 
ues obtained with the Hubble data for different 8 planes. These plots 
follows the same same order used in the two-dimensional C.L. of Fig.l 



Bounds for cosmic parameters 


9c 


Zt 


q(z = 0) 


Ueff(z = 0) 




Q z (z ~ 1100) 


I 


0.75 


-0.65 


-0.77 


0.02 


1.6 x 10~ n 


II 


0.66 


-0.55 


-0.70 


nn 9 


2.6 x 10~ 10 


III 


0.63 


-0.54 


-0.69 


0.01 


4.5 x HT 11 


IV 


0.75 


-0.62 


-0.74 


0.01 


4.3 x HT 11 


V 


0.68 


-0.60 


-0.73 


0.02 


2.4 x HT 11 


VI 


0.75 


-0.62 


-0.74 


0.01 


4.9 x 10~ n 



TABLE II: Derived bounds for cosmic parameters using the best 
fits value of 2D C.L. obtained in Table. ^ by varying two cosmolog- 
ical parameters in six different cases: I— (Ho^xo) — (74.32,0.77), 
//- (H a ,n m0 ) = (72.28,0.29), ///- (Ho, 7a) = (71.74, 10" 4 ), IV- 
(S! l0 ,!i m o) = (0.74,0.23), V- (fi^.To) = (0.74, 10~ 4 ), and VI- 
(ls,i1 m a) = (10~ 4 ,0.24). 



For the sake of completeness, we also report bounds for 
other cosmological relevant parameters [see Table (jlT|) ] . 
such as the fraction of radiation il r (z = 0), the effec- 
tive equation of state at z = (u; e //o = 7e//o — l)j the 
decelerating parameter at the present time qo, and the 
transition redshift (z t ); all these quantities are derived 
using the six best-fit values reported in Table fll]). We 
find that the z t is of the order unity varying over the in- 
terval [0.63, 0.75 ] ;su ch values are close to z t = 0.69+o;?° 
reported in [2J], [25| quite recently. Moreover, taking into 



account a \ 2 statistical analysis made in the (ajQ m ,z t ) 
plane based on the supernova sample (Union2), it has 
been shown that at 2a C.L. the transition redshift varies 
from 0.60 to 1.18 [26[. The behavior of a decelerating 
parameter with redshift is shown in Fig. 3, in particu- 
lar, its present-day value varies as —0.62 < q < —0.54 
for the six cases mentioned in Table II; all these values 
are in perfectly agreement with the one reported by the 
WMAP-7 project In Fig. © we plot the effective 
equation of state as a function of redshit for the best- 
fit value shown in Table (jTTJ) . In general, we find that 
— 1 < uj e ff < whereas their present-day values cover 
the range [—0.74, —0.69], so this does not exhibit a quin- 
tom phase [2l[. Regarding the behavior of density pa- 
rameters Q, x , f2 m , and f2 r , we find that nearly close to 
z = 0, the dark energy is the main agent that speeds up 
the Universe, far away from z — 1 the Universe is dom- 
inated by the dark matter and at very early times the 
radiation component enter in the action, controlling the 
entire dynamic of the Universe around z ~ 10 3 [cf. Fig. 
(J4j> ] . As it was expected the fraction of radiation at the 
present moment is negligible; thus, its value varies over 
the range 10~ 8 < Q r0 < 10~ 3 [sec Tabic (0)]. 

Now, we seek for another kind of constraint that 
comes form the physics at early times because this can 
be considered as a complementary tool for testing our 
model. As is well known, the fraction of dark en- 
ergy at the recombination epoch should fulfill the bound 
fl x (z ~ 1100) < 0.1 in order for the dark energy model 
be consistent with the big bang nucleosynthesis (BBN) 
data. Some signals could arise from the early dark en- 
ergy (EDE) models, uncovering the nature of DE as well 
as their properties to high redshift, giving an invaluable 
guide to the physics behind the recent speed up of the 
Universe p2l |. Then, the current and future data for con- 
straining the amount of EDE was examined, and the 
cosmological data analyzed has led to an upper bound 
of Q x (z ~ 1100) < 0.043 with 95% confidence level in 
the case of relativistic EDE, while a quintessence type of 
EDE has given n x (z ~ 1100) < 0.024 although the EDE 
component is not preferred, it is also not excluded from 
the current data |22||. Another forecast for the bounds 
of the EDE are obtained with the Planck and CMBPol 
experiments [1^; thus, assuming a Q x (a ~ 10~ 3 ) ~ 0.03 
for studying the stability of this value, it found that la 
error coming from the Planck experiment is (jP lanck ~ 
0.004 whereas the CMBPol improved this bound by a 
factor of 4 

Taking into account the best-fit values reported in 
Table (p}, we find that at early times the dark energy 
changes rapidly with the redshift z over the interval 
[10 3 , 10 4 ] [see Fig. g]) for more detail]. Indeed, Table (jn| 
shows that it varies as 2.4 x 10~ n < fl x < 2.6 x lO" 10 
around z ~ 1100. Such findings point out that the model 
constructed here not only fulfills the severe bound of 
£l x (z ~ 1100) < 0.1 but it is also consistent with the 
future constraints achievable by Planck and CMBPol ex- 
periments [2^] as well, corroborating that the value of the 
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cosmological parameters selected before, through the sta- 
tistical analysis made with Hubble data, are consistent 
with BBN constraints. 

It is interesting to compare the former analysis, where 
radiation was considered as a free evolving component 
which is decoupled from the dark sector in relation with 
other cases, where a radiation component interacts with 
the dark sector in order to provide a quantitative anal- 
ysis of the role played by the interaction. At present, 
dark energy dominates the whole dynamics of the Uni- 
verse and there is practically an obvious decoupling with 
radiation. However, from a theoretical point of view, 
it is reasonable to expect that dark components can in- 
teract with other fluids of the Universe substantially in 
the very beginning of its evolution due to a process oc- 
curred in the early Universe. For instance, dark energy 
interacting with neutrinos was investigated in (27| . The 
framework of many interacting components could pro- 
vide a more natural arena for studying the stringent 
bounds of dark energy at the recombination epoch. There 
could be a signal in favor of having dark matter exchang- 
ing energy with dark energ y w hile radiation is treated 
as a decoupled component [Til ]. jTEj or the case where 
dark matter, dark energy, and radiation exchange en- 
ergy More precisely, when the Universe is filled with 
an interacting dark sector plus a decoupled radiation 
term, it was found that Q x (z ~ 1100) = 0.01 [3 or 
VL x {z ~ 1100) = 10~ 8 [HI but if radiation is coupled to 
the dark sector, the amount of dark energy is drastically 
reduced, giving Q x (z ~ 1100) ~ ©(lO^ 10 ) so the ad- 
dition of an interacting radiation term is important for 
reducing the amount of early dark energy in 2 or 8 orders 
of magnitude. 



III. CONCLUSION 

We have taken under study the dynamical behavior of 
a cosmological scenario in which the Universe contains 
three interacting components, namely, dark energy, dark 
matter, and radiation within the framework of the usual 
spatially flat FRW spacetime. We have worked within 
the case of a transversal linear interaction because it 
gives a unique preferred direction in the plane constraint 
Yli Qi = ; then, for such a case we have obtained the 
partial energy densities (p Xl p m , p r ) in terms of the total 
energy density, and its derivatives up to second order. 
Additionally, we have imposed the existence of a power 
law solution and investigated its stability finding that the 
system is asymptotically stable in the sense of Lyapunov 
for a > and (3 > [see Eqs. (HH|) - ] ■ Using a linear 
transversal interaction that depends on the total density 
and its derivatives up to third order only, we have found 
that the interaction between the three components makes 
the total energy density exhibit three different behaviors 
with the scale factors a~ 3ls , a _37 + , and a _37 ~ , altering 
the usual corresponding traits to have no interaction. 

On the observational side, we have examined the pre- 



vious model by constraining the cosmological parameters 
with the Hubble data and the well-known bounds for dark 
energy at the recombination era. In the case of 2D C.L., 
we have made six statistical constraints with the Hub- 
ble function [see Fig. ([3]) and Table (U)]. We have found 
that j s varies from 10~ 4 to 10 -3 , so these values clearly 
fulfill the constraint j s < 2/3 for getting an accelerated 
phase of the Universe at late times. Using the priors 
(#0,7*) = (74.20km s" 1 Mpc" 1 , 10" 3 ), the best-fit val- 
ues for the present-day density parameters are given by 
(n x0 ,fl m o) = (0.74,0.23) along with a xl.o.f. = °- 779 < 
1. We have obtained that the estimated values of f2 m o 
vary from 0.24 to 0.29, without showing a significant dif- 
ference with the standard ones [2(|. Besides, it turned 
out that H G [71.28, 74. 32]km s" 1 Mpc" 1 so the latter 
values are met within la C.L. reported by Ricss et al 
flTjj . Regarding the derived cosmological parameters, for 
instance, the transition z t turned to be of the order unity, 
varying over the interval [0.63,0.75]; such values are in 
agreement with z t = 0.69tQ ^3 reported in (24|-[25[ , and 
meet within the 2er C.L obtained with the supernovae 
(Union 2) data in (26[. Besides, the decelerating param- 
eter q(z — 0) € [—0.62,0.54] and the equation of state 
ui e f f {z = 0) G [-0.74, -0.69]; indeed, -1 < u) eff < [sec 
Fig. (Q|], while the fraction of radiation at the present 
momment f2 r o varies in the interval [10 -3 , 10~ 9 ] for the 
six cases mentioned in Table ([II)). Further, the dark en- 
ergy amount 0, x (z) governs the dynamic of the Universe 
near z = 0, whereas far away from z = 1 the Universe is 
dominated by the fraction of dark matter f2 m (z) and at 
very early times the fraction of radiation Q r (z) controls 
the entire dynamic of the Universe around z ~ 10 3 [cf. 
Fig. ®]. 

Finally, taking into account the best-fit values reported 
in Table (P), we find that at early times the dark en- 
ergy changes rapidly with the redshift z over the inter- 
val [10 3 ,10 4 ] [cf. Fig. (@}]; indeed, we have obtained a 
2.4 x 10" 11 < fl x < 2.6 x 10~ 10 around z ~ 1100 [ see Ta- 
ble (JTTJ)] . The latter results indicate that the model con- 
structed fulfills the severe bound of Vt x {z ~ 1100) < 0.1 
and is consistent with the future constraints achievable 
by Planck and CMBPol experiments [23| as well, point- 
ing that the cosmological parameters selected before are 
coherent with BBN constraints also. In a future inves- 
tigation, we will perform a full study of the interaction 
vector proportional to orthonormal projection e Q that is 
related to the other direction defined in the constraint 
plane £\ Qi = 0. 
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